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Using a Navier- Stokes isotropic turbulent field numerically simulated in a box with a discretiza- 
tion of 1024^ (Biferale L. et al. Physics of Fluids, 17(2), 021701/1-4 (2005)), we show that the 
probability of having a stretching-tilting larger than twice the local enstrophy is negligible. By 
using an anisotropic kind of filter in the Fourier space, where wavenumbers that have at least one 
component below a threshold or inside a range are removed, we analyze these survival statistics 
when the large, the small inert ial or the small inertial and dissipation scales are filtered out. It can 
be observed that, in the unfiltered isotropic field, the probability of the ratio (|u? • VU|/|u?|^) being 
higher than a given threshold is higher than in the fields where the large scales were filtered out. At 
the same time, it is lower than in the fields were the small inertial and dissipation range of scales 
is filtered out. This is basically due to the suppression of compact structures in the ranges that 
have been filtered in different ways. The partial removal of the background of filaments and sheets 
does not have a first order effect on these statistics. These results are discussed in the light of a 
hypothesized relation between vortical filaments, sheets and blobs in physical space and in Fourier 
space. The study in fact can be viewed as a kind of test for this idea and tries to highlight its limits. 
We conclude that a qualitative relation in physical space and in Fourier space can be supposed to 
exist for blobs only. That is for the near isotropic structures which are sufficiently described by a 
single spatial scale and do not suffer from the disambiguation problem as filaments and sheets do. 

Information is also given on the filtering effect on statistics concerning the inclination of the strain 
rate tensor eigenvectors with respect to vorticity. In all filtered ranges, eigenvector 2 reduces its 
alignment, while eigenvector 3 reduces its misalignment. All filters increase the gap between the 
most extensional eigenvalue < Ai > and the intermediate one < A2 > and the gap between this 
last < A2 > and the contractile eigenvalue < A3 >. When the large scales are missing, eigenvalue 
modulus 1 and 3 become nearly equal, similar to the modulus of the related components of the 
enstrophy production. 

PACS numbers: 47.27.Ak, 47.27.eb, 47.27.Gs, 47.32.C- 



I. INTRODUCTION 

The formation of spatial and temporal internal scales 
can in part be associated to the stretching and tilting of 
vortical structures. Many aspects of the behavior of tur- 
bulent fields have been associated to this phenomenon: 
the onset of instability, vorticity intensification or damp- 
ing, or the three-dimensionalization of the flow field [WS] . 
In the standard picture of turbulence, the energy cascade 
to smaller scales is interpreted in terms of the stretching 
of vortices due to the interaction with similar eddy size 
(see for example [4 ). A number of statistical details on 
the stretching phenomenon and the closely related en- 
strophy production can be found in the monography by 
Tsinober (2001, see in particular Chapter 6, [5]). 

Although the important physical role of these inertial 
phenomena is recognized, the literature does not often 
include statistical information on quantities such as the 
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magnitude or the components of u) • VU. For instance, 
in a letter to Nature (2003) dedicated to the measure- 
ments of intense rotation and dissipation in turbulent 
flows, Zeff et al. [6 observe that the understanding of the 
temporal interactions between stretching and vorticity is 
crucial to the science of extreme events in turbulence. 
However, the statistics presented there concern dissipa- 
tion and enstrophy and not directly stretching. The lit- 
erature more often includes statistical information con- 
cerning other gradient quantities such as the strain rate 
or the rate-of-rotation tensors, and, in particular, their 
fundamental constituents: the longitudinal or transverse 
velocity derivatives. Over the last 20 years, statistics on 
the skewness and flatness factors of the velocity deriva- 
tive have been considered in a number of laboratory and 
numerical studies that show how these quantities increase 
monotonically with the Reynolds number, see e.g. [7 and 
the review by Sreenivasan and Antonia (1997) [8|. 

In the case of turbulent wall flows, laboratory measure- 
ments of both the mean and the r.m.s. of fluctuations 
of the stretching components across the two-dimensional 
boundary layer have been reported by Andreapoulos and 
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Honkan (2001) [9]. In this study, the normahzed r.m.s 
values of the stretching components are very significant 
throughout the boundary layer and reach values that are 
one order of magnitude larger than the mean span-wise 
component (the only significant mean component, how- 
ever and only in the near wall region). The values ob- 
served for the r.m.s. of the stretching range from 0.04, 
close to the wall, to about 0.004 in the outer part. 

In a study concerning the structure and dynamics of 
vorticity and rate of strain in incompressible homoge- 
neous turbulence, Nomura and Post (1998), [10 demon- 
strate the significance of both local dynamics (infiuence 
of local vorticity) and spatial structure (infiuence through 
non-local pressure Hessian) in the interaction of the vor- 
ticity and strain rate tensor. The behaviour of high- 
amplitude rotation-dominated events cannot be solely 
represented by local dynamics due to the formation of 
distinct spatial structure. Instead, high-amplitude strain 
dominated regions are generated predominantly by lo- 
cal dynamics. The associated structure is less organized 
and more discontinous than the one associated with ro- 
tation dominated events. They conclude that non-local 
effects are significant in the dynamics of small scale mo- 
tion. This should be considered in the interpretation of 
single-point statistics. Characterizations of small-scale 
turbulence should consider not only the typical structures 
there present but also typical structure interactions. In 
this context these authors offer the radial distribution of 
the magnitude of the strain rate tensor normalized on 
the enstrophy. In this paper the maximum value of this 
magnitude is found close to 0.2. 

Laboratory statistical information on the stretching of 
field lines can be found in [H . Here, probability density 
functions of the logarithm of the local stretching in N 
cycles were obtained for several two-dimensional time- 
periodic confined fiows exhibiting chaotic advection. The 
stretching fields were observed to be highly correlated in 
space when N is large, and the probability distributions 
were observed to be similar for different fiows. 

However, a few examples in literature can also be 
cited regarding direct results for stretching-tilting statis- 
tics. For instance, recently experimental and numeri- 
cal confirmation has been found of the predominance 
of three dimensional turbulent vortex stretching in the 
positive net enstrophy production. These aspects have 
been extensively considered in Tsinober (2000) [12] and 
in the 2001 monography [5j, where a number of statisti- 
cal geometrical details concerning the vortex alignement, 
compression, tilting, and folding are outlined. Through 
two papers, Constantin, Procaccia and Segel (1995) [13 , 
Galanti, Procaccia and Segel (1996) |14 consider the 
stretching and its relationships with the amplification of 
vorticity and the straightening of the vortex lines. They 
show that the same stretching that amplifies the vorticity 
also tends to straighten out the vortex lines. They also 
show that in well-aligned vortex tubes, the self-stretching 
rate of the vorticity is proportional to the ratio of the 
vorticity and the radius of curvature. In this context, [M] 



gives statistics on the stretching and vortex line curva- 
ture. Numerically this is seen as the appearance of high 
correlations between the stretching and the straightness 
of the vortex lines. Regarding to this issue, an impor- 
tant universal feature of fully developed turbulent fiows 
is the preferential alignment of vorticity along the eigen 
direction of the intermediate eigenvalue of the strain-rate 
tensor. A number of works both experimental and nu- 
merical studies on this result are available (Tsinober, Kit 
and Dracos JFM (1992), Kholmyansky, Tsinober and 
S. Yorish PoF (2001), [l5], Gulitski et al. JFM (2007 
a,b,c), [HH^ and Chevillard et al. (2008), [ig). It 
should be noticed, however, that in the case of nonlo- 
cal strain rate, Hamlington, Schumacher and Dahm [20] . 
have observed a direct assessment of vorticity alignment 
with the most extensional eigenvector by using data from 
highly resolved direct numerical simulations. 

In the present study, for the case of isotropic turbu- 
lence {Rex = 280, [21 ), we consider statistics related 
to the intensity of the stretching term in the equation 
for vorticity. If we consider the general instantaneous 
local intrinsic anisotropy of turbulent fields, looking at 
stretched structures as filaments and sheets, we would 
like to be able to disentangle them to follow and un- 
derstand better their evolution and detailed dynamics. 
Isotropic filtering is unable to carry out this job. 

We have conceived a probe function, the ratio between 
the magnitude of the vortex stretching and the enstro- 
phy, to empirically and statistically measure the local 
activity of the stretching phenomenon (see section II). 
In addition, we propose an alternative to the commonly 
used isotropic filter: the cross filter. This is a new em- 
pirical, and at the moment limited, attempt to intro- 
duce an anisotropic filtering. In section HI, we analyze 
the survival function of the normalized stretching by us- 
ing the cross filter acting directly on the velocity Fourier 
space. We do this in the hope of qualitatively highlight- 
ing aspects related to the role of the three-dimensional 
structures known as blobs, sheets and filaments and their 
hypothetical Fourier counterparts. This study can be 
viewed as a kind of test for this idea and tries to high- 
light its limits. Concluding remarks are made in section 
IV. 



II. THE NORMALIZED 
STRETCHING-TILTING FUNCTION 

With reference to the phenomena described by the in- 
ertial nonlinear nonconvective part of the vorticity trans- 
port equation, let us introduce a local measure of the 
process of three-dimensional inner scales formation 



/(x,^) = ^(x,t) = ^(x,t). (1) 



where U is the velocity field, Sij is the strain rate ten- 
sor, and uj is the vorticity vector. The numerator, the so 
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called stretching-tilting term of the vorticity equation, is 
zero in two-dimensional flows. In 3 D fields, it is com- 
monly believed to be responsible for the transfer of the 
kinetic energy from larger to smaller scales (positive or 
extensional stretching) and viceversa (negative or com- 
pressional stretching). According to definition (1), / de- 
pends on the local instantaneous velocity and vorticity 
fields. In this study, we leave aside the peculiarity asso- 
ciated to the convective forcing and focus on the action 
of the fiuctuation field only. For simplicity, we consider 
here the fiuctuation of an homogeneous isotropic turbu- 
lent field (^J). Since the stretching term plays an im- 
portant role in the enstrophy production, in the previ- 
ous definition the normalization by \u)\^ was adopted. It 
should be recalled that the square of the vorticity mag- 
nitude is the only invariant of the rate of rotation tensor 
which is non zero and is also the square of the Frobenius 
norm, an invariant norm of the rate of strain tensor. For 
this reason, we considered the enstrophy a good candi- 
date as reference quantity for the product u) • VU. In 
fact, as it can be seen below, the survival probability dis- 
tribution function of / is very small for values larger than 
0(1). 

Function / was evaluated over a fully resolved homo- 
geneous isotropic incompressible steady in the mean tur- 
bulence in order to look for the typical range of values 
of /(x) and to relate them to the behavior of the various 
turbulence scales present in an isotropic field. 

The dataset consists of 1024^ resolution grid point Di- 
rect Numerical Simulation (DNS) of an isotropic Navier- 
Stokes forced field at Reynolds Rex = 280 [21]. Each 
instant in the simulation is statistically equivalent, and 
provides a statistical set of a little more than 10^ ele- 
ments. We considered the statistics that were obtained 
averaging over the full domain in one instant. The field 
has been slightly modified in order to filter out instanta- 
neous effects of the forcing, in other words, a turbulent 
kinetic energy inhomogeneity of about 20% (in the spa- 
tial coordinate system). As this bias was generated by 
the energy supply at the large scale range, the two largest 
scales have been filtered out. The resolved part of the en- 
ergy spectrum extends up to /c ~ 330. The inertial range 
extends from /c ~ 10 to A: ~ 70, see the compensated 
version of the 3D spectrum in figure [l] The higher wave- 
numbers, which are affected by the aliasing error, are not 
shown. 

We focus now on a few statistical properties of . 
The pdf of the components of this vector (which are sta- 
tistically equivalent, since the field is isotropic) is shown 
in figure [2] Symmetry with the vertical axis is expected 
because of isotropy; the skewness is in fact approximately 
10~^, which is not meaningfully far from zero. However, 
the distribution cannot be approximated with a Gaussian 
function. In fact, the actual kurtosis is approximately 55, 
which is very far from the Gaussian value of 3. 

The range of values attained by /(x) is wide. Val- 
ues as high as a few hundreds were observed at a sparse 
spatial points. In order to read the typical values of 
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FIG. 1: Compensated 3D energy spectrum of the isotropic 
field here considered, Re\ — 280 
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FIG. 2: Probability density function of one component of 
the vector ^^2^ in an isotropic velocity field with Re\ — 
280. Comparison with the Gaussian model. The Skewness 
is negligible, about 10~^. The Kurtosis however is very high 
and reaches a value of about 55. 



/(x), we study its survival function. By denoting F{s) = 
P{f{^) ^ the cumulative distribution function (cdf) of 
/(x), the survival function is defined as the complement 
to 1 of the cdf, 

S{s) = P{f{x) >s) = l- F{s). (2) 

For each threshold 5, S{s) describes the probability that 
/(x) takes greater values than s. 

It has been found that, when /(x) is evaluated on a well 
resolved isotropic turbulent field, the probability that 
/(x) > 2 is almost zero, see figure [s] Thus, /(x) = 2 can 
be considered the maximum statistical value that /(x) 
can reach when the turbulence is fully developed. 
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FIG. 3: Survival probability of the normalized stretching- 
tilting function in a fully resolved isotropic 3D turbulent field 
(^(/(x) > s), Rex = 280. Unfiltered velocity field. The 
dashed vertical line indicates the value of / where the proba- 
bility density function is maximum. 



III. PROPERTIES OF THE SURVIVAL 
FUNCTION OF THE NORMALIZED 
STRETCHING-TILTING TERM: ANALYSIS ON 
THE ANISOTROPICALLY FILTERED FIELD 

By means of suitable convolutions, the application of 
filters to the velocity field allows the behavior of the func- 
tion f{x) to be studied in relation to the different turbu- 
lence scale ranges. This analysis is carried out using two 
spectral filters, a high pass and a band-stop filter. To 
focus in an empirical way on the three principal kind of 
geometrical structures observed in turbulence, filaments, 
sheets and blobs, we use here a highly anisotropic kind 
of filter, which is less traditional than the axisymmetric- 
type filter. Of course, given the inadequacy of the spec- 
tral representation to account for the complex three- 
dimensional geometry of the turbulent structures com- 
monly seen through visualization tools, the approach we 
use here is not rigorous and should be considered no more 
than propaedeutical. 

The first filter is a sort of high-pass filter, which we 
refer to as cross filter and which allows the contribution 
of the structures that are characterized by at least one 
large dimension to be removed. From the Fourier point 
of view, this means that the structures whose wave- vector 
has at least one small component are filtered out. One 
can here think about elongated structures as filaments 
and sheets or very large globular structures. In figure |4j 
a graphical scheme of the filtering in the wave number 
plane /ci,/c2 can be seen. The first filter we consider is 
here represented in blue and is a kind of high-pass filter 
which affects all wave- numbers that, along any possible 
direction, have at least one component under a certain 
threshold. Given the threshold /cm/at, the filter reduces 
the contribution of the modes with wave number compo- 
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FIG. 4: Schema of the anisotropic filter here named as CROSS 
filter. Blue region: high-pass filter, the wave-numbers under 
a certain threshold are partially removed, see eq. (4), Red 
region: band-stop filter, the wave-numbers inside a range are 
cut, see eqs. (5-7). 
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ki < kuiN or /c2 < kuiN or < kuiN- 

The representation of this high-pass filter, g^p^ can be 
given by a function of the kind 

9hp{K) = Y[(l){ki;kMiN), (t){ki, kuin) = ^ ^ ^-(j^^-Um!^^ 

Since function g^p filters any wavenumber that has at 
least one component lower than the threshold kuiNi it 
reduces the kinetic energy of the filamentous (one compo- 
nent lower than kMiN)^ layered (two components lower 
than kMm) and blobby (three components lower than 
kuin) structures. This filter is efficient in reducing the 
integral scale of the turbulence [22 . 

By varying the value of the threshold, /cm/at, it is 
possible to consider different scale ranges. The ranges 
— 10, — 20, — 40 are compared in figure [5j The first 
filtering affects the energy-containing range, while the 
other two also include a part of the inert ial range, which 
is visible in figure [l] 

The plots in figure 5 have coherent behavior. The sur- 
vival function S for the — 10 filtering is slightly be- 
low the values of the distribution of the unfiltered turbu- 
lence. This trend is confirmed by the other two filterings, 
and the reduction grows as the threshold kuiN increases. 
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FIG. 5: Survival probability of the normalized stretching- 
tilting function in a high pass filtered isotropic turbulent field. 
CROSS filter, see in fig. 4 the blue region. 



The high-pass filter has the effect of decreasing the sta- 
tistical values taken by f{x) in the domain. The wider 
the filtered range, the higher the effect on /. 

It is possible to say that when we reduce the weight of 
the large-scale structures (layers, filaments or blobs), the 
local stretching-tilting intensity decreases with respect to 
the vorticity magnitude. On average, the values of f{x) 
go down. The wider the range affected, the lower the 
probability value becomes. This suggests that the large 
scales contribute more to the stretching-tilting (the nu- 
merator of /) than to the magnitude of the vorticity fluc- 
tuation (the denominator of /). It should be noted that 
this trend is consistent with the results in [20^ . This con- 
sistency also includes results relevant to the behaviour of 
the stretching fluctuation and of the vorticity magnitude 
in boundary layer turbulence, see figures 6 and 9 in [9]. 
For the wider range — 40, a decrease of 30% in the cu- 
mulative probability is observed for a stretching-tilting of 
about one half of the local vorticity. The decrease goes 
up to 80% when statistically the stretching-tilting has 
the same magnitude of the vorticity, that is / is close to 
1, see figure [5] 

Let us now consider the behavior of f{x) when the 
inertial and dissipative ranges are affected by the filter- 
ing, namely a band-stop filtering. In this case, the band 
width can be extended to obtain a low pass filtering. 

This filter can be obtained by reducing the contribu- 
tion of a variable band (see figure lij part in red) 



^MiN < h < kuAX or kuiN < ^2 < kuAX or 
kuiN < ks < kuAX- 



FIG. 6: Probability of the normalized stretching-tilting func- 
tion in a band pass filtered isotropic turbulent field of being 
higher than a threshold s. Control function: survival function 
1 — F(x), band-stop filtering in various portion of the inertial 
range and in the dissipative range. 

This yields the filter function 

9hs{K) \\4>{h]kMiN,kMAx), (4) 

(t){ki]kMIN,kMAx) = - (l>{ki]kMIN)\^ (l>{ki]kMAx) 

The effects of the application of this band-stop filter 
on the probability P{f{'x.) > s) are shown in figure [oj 

Let us now consider the inertial range in an extended 
way, which includes the — | range plus all the scales 
which are not yet highly dissipative. The different bands 
are 10 — 40, large scale inertial filtering, 40 — 70 interme- 
diate scale inertial filtering, 70 — 100 small scale inertial 
filtering, 100 — 130 near dissipative, 30—150 intermediate- 
inertial/small scale filtering, 150 — 330 dissipative scale 
filtering. Once again all the filtered ranges induce the 
same effects: for 5 < 1/2, a slight increase in the survival 
probability. For small values of 5, the most effective fil- 
tering (i.e. the ones which produce the highest increase) 
is the 10 — 40, while, for higher statistically relevant val- 
ues, 0.5 < s < 2, the most effective result is obtained 
filtering over the whole inertial range, 30 < A: < 150. In 
this case, an increase of about 60% is observed for 5 = 1 
and of about 80% for s = 1.5 

This highlights the fact that the structures of the in- 
ertial range contribute more to the intensity of the vor- 
ticity field than to stretching and tilting. The general 
trend is almost inverted with respect to the case of the 
high pass filtered turbulence (compare the 0-40 and 10-40 
results in figures 4 and 5, respectively) and this can be 
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confirmed, with sligiit differ ences, as long as we enlarge 
the amplitude of the filtering band to get closer to the 
dissipative range. Finally, moving toward the dissipative 
range (150 < k < 330), the band-stop filter becomes a 
sort of low-pass filter. By filtering these wave numbers, 
the obtained effect is minimum, although we have re- 
moved the contribution of more or less the highest 200 
wave-numbers (see figure [6|. 

To see the effect of the filters on the vortical structures, 
the vorticity magnitude has been visualized in two ways. 
The first is a volume rendering of the surfaces where one 
of the vorticity components is close to the root mean 
square value, see figure [7[ The instantaneous field we 
are considering is homogeneous and isotropic, so all the 
components are statistically alike and it suffices to ob- 
serve one component only. In panel (a) the unfiltered 
field is visualized and a complex picture made of an elon- 
gated, thick, sleeve-like structures which are enfolded and 
twisted can be seen. In panel (b) the structures with at 
least one wavenumber component below 20 are smoothed 
out, see equation ([3|. One here sees a more sparse dis- 
tribution of mostly elongated and nearly flat structures 
with a much shorter length with respect to panel (a) . The 
mutual folding and twisting seems reduced. This image 
is related to the survival distribution in figure 5, where 
a depression of stretching-tilting over the vorticity mag- 
nitude is reported for ratio values above 0.3. Panel (c) 
shows the band-stop filtered field, where wavenumbers in 
between 30 and 150 are smoothed out. Basically, most 
of the inertial and larger dissipative structures are re- 
moved. Here, it is interesting to observe that the surface 
is much more corrugated that in panel (a), which leaves a 
definitive view of the small scale above wavenumber 150. 
Some of these structures are elongated, others are glob- 
ular. The image does not discourage the idea that the 
large and the very small scales directly interact. Finally, 
in panel (d) one sees the structures that have at least 
one wavenumber component in the range 0-30. Here, 
large unruffled structures which are mainly globular can 
be seen. The images in panel (c) and (d) represent in- 
stances of the statistical situation described in figure 6, 
where the survival ratio of the stretching and vorticity 
intensities is enhanced with respect to the natural sit- 
uation. Thus, it seems that the partial absence of the 
inertial range amplifies the stretching-tilting process. It 
should also be noted that the structure spatial density 
is distributed in almost the same way in all the panels, 
though the density levels are different. 

We have tried to visualize the filtering effect also by 
means of contour plots and pseudo color imaging of the 
vorticity magnitude in a flat section of the 1024^ fields see 
figure [S] Here, the lines in the images in the left column 
are the contours of iso-surfaces of the vorticity magni- 
tude. Starting from the top, one sees the unfiltered field, 
the high pass filtered field (the wavenumbers above 20 are 
kept), the 30 — 150 band-stop, and the low pass filtered 
field (the range 30 — infinity is removed). The contour 
plot technique is very popular, but, apart from clearly 



showing the reduction of the turbulence scale size when 
the large scales are missing, it does not give much infor- 
mation. Essentially, structure contours appear to have 
the peanut shape which typically hosts vortex dipoles. 

In the enlarged views in the central and right columns, 
the pseudocolor plots are richer in information. Where 
the larger scales are removed, panels (e) and (f), the sur- 
vival probability of intense stretching is reduced. The 
range of variation for the vorticity magnitude is very 
large, the root mean square value is about seven times 
smaller than the maximum value. When only the large 
scales are left, panels (m) and (n), the rms value is about 
one third of the maximum value and the stretching is en- 
hanced. When the inertial scales are removed, the large 
scale appears to be wrapped by the small scales. The big 
dipoles are surrounded by thin wavy-like sheets and the 
small scales are attached to the large ones. In this situa- 
tion it is not possible to neglect their direct interaction. 
In regions where a large scale is missing, the small scales 
are also missing. Viceversa, a low level of direct inter- 
action between the largest and the smallest scales would 
have been confirmed, if small scales had been sparsely dis- 
tributed in regions where the large scales are not present 
and would not have surrounded the large structures in 
regions where the latter are present. 

To complement the understanding of the visualization 
in relation to the anisotropic filtering here used, we have 
considered the alignment between the eigenvectors of the 
strain rate tensor and the direction of the vorticity. Func- 
tion / in fact, see equation (1), can also be written as 

/(x) = j^\[sf{ei * Gu;)^]"^/^ where l^l is the magnitude 
of the strain rate tensor, Si are the eigenvalues of the 
strain rate tensor Sij normalized hy \S\^ Si = Xi/\S\, and 
i^i * ^uj) describe the alignments between the eigenvec- 
tors of Sij, denoted e^, and the direction of the vortic- 
ity . Figure [9] shows the probability density function 
of these alignments in the reference unfiltered field and 
in two filtered cases described in this work: the high 
pass filter where the smallest 20 wavenumbers are re- 
moved, see equation (3) and the band-stop filter where 
wavenumbers 30 — 150 are removed, see equation (4). 
We first observe that the standard trend of alignments is 
not fully spoiled by the filtering. In both filtering cases, 
eigenvector 2 reduces its alignment, while eigenvector 3 
reduces its misalignment. Conversely, eigenvector 1 in- 
stead shows a different behavior. In the band-stop fil- 
tering case (large scale dominate) eigenvector 1 slightly 
increases the alignment. In the high-pass filtering, eigen- 
vector 1 reduces the alignment that becomes statistically 
equal to that of the eigenvector 3. This is confirmed, 
see Table 1, by considering the ratio among the field av- 
eraged strain rate tensor eigenvalues and related com- 
ponent of the enstrophy production, < > / < atot 
where < • > is the average over all the computational do- 
main (1024^ grid point box), = u'^ Xicos'^ {e^ , ei) and 
o'tot =< ^'^ Xicos'^ {e^ ^ Ci) >. All filters increase the 
gap between the eigenvalue < Ai > and < A2 > and the 
gap between < A2 > and < A3 >. However, when the 
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FIG. 7: Three dimensional visualization of the surfaces where one vorticity component has the value 17. 5sec~^. The root 
mean square value of the vorticity magnitude in the field is |c(;|rms = 30.3sec~^. As reference, we consider a water field with 
Rex = 284, Urms = 0.09m/ sec, v — 0.96 * \0~^m^ j sec (water viscosity at 23 °C), Taylor micro-scale A = 3mm and integral 
scale / = 56mm. Panel (a): unfiltered field; panel(b): the wave number range 0-20 is filtered out by using the high-pass 
cross filter, panel (c): the wave number range 30-150 is filtered out by using the band-stop cross filter; panel (d): the wave 
number range 30-infinity is filtered out by using the low-pass cross filter, i.e. by letting kuiN oo, see figure [4] The 
visualization shows a 256^ portion of the numerical field simulated on a 1024^ point grid. The field is visualized by means of 
Visit (https:/ /wci. llnl.gov/codes/visit/). 



large scales are missing, eigenvalues 1 and 3 are very close 
in modulus. The same happens to the modulus of their 
related enstrophy production components. When the in- 
ertial scales and part, or the entire, dissipative range are 
removed the mutual relation among the eigenvalues and 
the modulus of the production components changes less 
with respect to the natural turbulence. 

In all the three cases the filtering reduces the average 
values of the largest and intermediate eigenvalue, < Ai > 
and < A2 >. 

At this point, let us consider the dual nature of the fil- 
aments and sheets, as regards their inclusion in the cat- 
egories of the small and large scales. A filament which 



is filtered out by the filter ghp because it has a small 
wave number (the axial wave number component), will 
have two large wave number components (the ones nor- 
mal to the filament axis). Due to these wave components, 
it will also be filtered out by the filter ^55. A similar 
situation also holds for the sheets. Thus filaments and 
sheets are always partially removed when a filtering, ei- 
ther anisotropic or isotropic, is applied. The situation 
changes with compact structures (the blob), which non 
ambiguously belong either to the large scale range or to 
the intermediate-small scale range. The different behav- 
ior shown in figures |5] and |6] is therefore mainly due to 
the blob contributions, and, since the variation in the cu- 



8 



Contour 
r^Var: MV_multi00_15 




PssLcbcofor 
Var:MV mu»tDO 15 
—-"436.1 " 




I- 17.22 
13.78 
10,33 
3.444 
O.DOO 



0.003555 
Max: 436.1 
Mln: 0.003535 




fteudscobr 
Var:MV muttDO 15 

^_-'3o.oa " 



Rseudocobr 
Var; MV multDO 15 
^H-'223.e ~ 



I Reinocotor 
I Van MV multDO 15 
^190,5 



I 



-0.003437 
■ Max: 190.5 
I Mln; 0.003437 



I Pseudocobf 
I Var:MV muttDO 15 
^-"96.41 



I 



-0,01794 
: 96.41 
■ l.lln: 0.01794 





fteudocobr 
Var: MV miitDO 15 



10.00 
Max: 223.8 
Mln; 0.01098 



fteudocobr 
Vaf:MV mmtDO 15 
^^-"30.00 ~ 



10.00 
Max: 190,5 
Mln: 0.003437 



fteudocobr 
Var: MV multDO 15 
^.00 - 



10,00 
M!.] . : 96.41 
Mh: 0.01794 



FIG. 8: Visualization of the vorticity magnitude in a section parallel to one face of the computational box. Reference isotropic 
turbulence: Rex = 280, Urms = 0.09m/ sec (rms of the velocity), v = 0.96 * 10~^ / sec (water viscosity at 23 ^C), Taylor 
micro-scale A = 3mm and integral scale / = 56mm, |cc;| rms = 30.3sec~^ (root-mean-square of the vorticity magnitude), 1024^ 
grid domain points. First row (a,b,c): unfiltered field. Second row (e,f,g): the wave number range 0-20 is filtered out by using 
the high-pass cross filter. Third row (g,h,i): the wave number range 30-150 is filtered out by using the band-stop cross filter. 
Fourth row (l,m,n) the wave number range 30-infinity is filtered out by using the low-pass cross filter. It is possible to see: 
in the first column (a,d,g,l) the vorticity magnitude countourplots of the entire 1024"^ grid domain. In the second and third 
column (b,e,h,m and c,f,i,n) Pseudocolor plots of a 256^ portion of the grid (black box in the previous column) are shown. In 
the third column the range of magnitude values in between 10 and 30 is visualized to show the details of the part of the field 
where the vorticity is below its rms value. The part above the rms is in fact visualized in the central column, where all the 
range of values is included. The Pseudocolor method maps the data values of a scalar variable to color. The plot then draws 
the colors onto the computational mesh. The field is visualized by means of Visit (https:/ /wci. llnl.gov/codes/visit/) 
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TABLE I: normalized strain rate tensor eigenvalues and normalized enstrophy production computed in 
the unfiltered and the filtered field 



Normalized strain rate tensor eigenvalues <Aj >/(<Af + 


+ A3 >) 


Principal 


Unfiltered 


Inertial 


Large and 


Large scales 
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index axes 


field 
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dissipative scales 


*** 


field 
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dissipative 


** 
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scales 
* 






L 1 oj 


I 


0.54 


0.6121 


0.57 


0.56 


0.47 


II 


0.12 


0.0013 


0.091 


0.090 


0.06 


III 


-0.66 


-0.6108 


-0.64 


-0.64 


-0.53 




Normalized 


enstrophy production <ai>/<atot> 




I 


1.14 


6.31 


2.67 


2.86 


1.06 


II 


0.66 


0.59 


0.49 


0.45 


0.51 


III 


-0.80 


-5.91 


-2.16 


-2.31 


-0.57 



* the smallest 20 wavenumbers are filtered out 

** wavenumbers in the inertial 30-150 range are filtered out 

*** all wavenumbers above 30 are filtered out 



i=3 unfiltered field — 




) I 1 1 1 1 

-1 -0.5 0.5 1 



a better way to select and remove anisotropic structures 
such as filaments and sheets is found, first order statisti- 
cal modification associated to their presence/absence will 
not be clearly seen. 

Lastly, it is interesting to observe that box filtering 
small scales modifies the stretching statistics to a great 
extent. A field filtered in such a way shows a finite 
probability of having more than twice the amount of 
stretching/tilting compared to the enstrophy [23, 24 . In 
the context of the Large Eddy Simulation methodology, 
where this kind of filtering is commonly used, it is pos- 
sible to deduce that, when a fluctuating field shows such 
a feature, the field is unresolved. As a consequence, it 
is possible to build a criterion that locates the regions 
of the field where the inclusion of a subgrid term in the 
governing equations is advisable (see also [25 J. 



FIG. 9: PDFs of the cosine of the angle between vorticity 
(u?), and the eigenvector, e^, of the rate of strain tensor (Sij). 
The red lines refer to unfiltered field; the green lines refer to 
the wave number range 0-20 filtered out by using the high- 
pass cross filter; the blue lines refer to the wave number range 
30-150 filtered out by using the band- stop cross filter. 



mulative distributions is opposite and almost of the same 
magnitude, it is possible to deduce that the partial re- 
moval of the background filaments and sheets, which is 
always done regardless of the filter typology (high pass, 
band pass, low pass,..), is not statistically relevant. In 
other words, it appears reasonable to conclude that until 



IV. CONCLUSIONS 

To summerize, we have collected a set of statisti- 
cal information about the stretching and tilting inten- 
sity of vortical structures normalized by the enstrophy, 
/(x) = ^^^^^ (x), in isotropic turbulence. A first re- 
sult is that there is a very small probability of having a 
larger stretching/tilting of intensity than the double of 
the square of the vorticity magnitude. Then, if compact 
structures (blobs) in the inertial range are filtered out, it 
can be seen that the probability of having higher / than 
a given threshold s increases by 20% at s = 0.5, and by 
60-70% at 5 = 1.0. If, on the other hand, larger blobs 
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are filtered, an opposite situation occurs. Tiie unfiltered 
field is thus a separatrix for the cumulative probability 
function. This behavior - high fluctuation vorticity mag- 
nitude low stretching, and viceversa - agrees with gen- 
eral aspects highlighted by a number of laboratory and 
numerical analyses, [TTl [20] also in near wall turbulent 
flow configurations [T. The present observations need to 
be associated to the non discriminating effect of filtering 
on filaments and sheets, which is due to their specific na- 
ture that cannot be reconciled inside either a category 
of small or large scales. It has also been shown that a 
high intermittency is associated to /, whose kurtosis is 
as high as 55. 

The probability density function of the alignments be- 
tween the strain rate eigenvectors and the vorticity is 
in part modified by the anisotropic filtering here investi- 
gated. In particular, we observe that, though the stan- 
dard trend of alignments is not fully spoiled, eigenvec- 
tor 2 reduces its alignment, while eigenvector 3 reduces 
its misalignment. Conversely eigenvector 1 shows a dif- 
ferent behaviour. In the band stop filtering case (large 
scale dominate) eigenvector 1 slightly increases the align- 
ment. In the high pass filtering case (inertial scales dom- 



inate), eigenvector 1 reduces the alignment that becomes 
statistically equal to that of the eigenvector 3. This is 
confirmed by considering the mutual ratio among the av- 
eraged strain rate eigenvalues and related components of 
the enstrophy production. Both filters increase the gap 
between the most extensional eigenvalue < Ai > and the 
intermediate one < A2 > and the gap between this last 
< A2 > and the contractile eigenvalue < A3 >. How- 
ever, when the large scales are missing, the modulus of 
eigenvalues 1 and 3 become nearly equal, similar to the 
modulus of the related components of the enstrophy pro- 
duction. 



V. ACKNOWLEDGMENTS 

We are grateful to Professor James J. Riley for fruitful 
discussions. We thank the referees for their constructive 
reviews and the many suggestions. This work was carried 
out in cooperation with the International Collaboration 
for Turbulence Research. 



[1] A.S. Monin and A.M. Yaglom Statistical Fluid Me- 
chanics. Vols 1-2, The MIT Press, Cambridge, Mas- 
sachusetts, and London, England., (1971, 1975). 

[2] H. Tennekes and J.L. Lumley A first course in turbulence, 
The MIT Press, Cambridge, Massachusetts, and London, 
England., (1972). 

[3] S.B. Pope, Turbulent flows, Cambridge University Press, 
(2000). 

[4] U. Frisch, Turbulence - The legacy of A.N. Kolmogorov, 
Cambridge University Press, (1995). 

[5] A. Tsinober, An Informal Introduction to Turbulence, 1th 
edition, Springer, (2001). 

[6] B.W. Zeff, D. Lanterman, R. McAllister, R. Roy, J. 
Kostelich and D.P Lathrop , Nature 421, 146-149 (2003). 

[7] A. Tsinober, E. Kit and T. Dracos, "Experimental in- 
vestigation of the field of velocity gradients in turbulent 
flows", J. Fluid Mech. 242, 169-192 (1992). 

[8] K. R. Sreenivasan and R. A. Antonia, Annu. Rev. Fluid 
Mech. 29, 435 (1997). 

[9] Y. Andreapoulos and A. Honkan, J. Fluid Mech., 439, 
131 (2001). 

[10] K.K. Nomura and G.K. Post, J. Fluid Mech. 377, 65-97 
(1998). 

[11] PE. Arratia and J.P Gollub, J. Stat. Phys., 121, Nos. 
516, 805 (2005). J. Fluid Mech., 439, 131 (2001). 

[12] A. Tsinober Turbulence structure and Vortex Dynamics, 
J.C.R. Hunt nad J.C. Vassilicos editors, Cambridge Uni- 
versity Press, 164 (2000). 

[13] P. Constantin, I. Procaccia and D. Segel, Phys. Rev. E. 
51, 3207 (1994). 



[14] B. Galanti, I. Procaccia and D. Segel, Phys. Rev. E. 54, 
5122 (1996). 

[15] M. Kholmyansky, A. Tsinober and S. Yorish, Physics of 

fluids, 13, 311 (2000). 
[16] G. Gulitski, M. Kholmyansky, W. Kinzelbach, B. Lthi, 

A. Tsinober and S. Yorish, J. Fluid Mech. 589, 57-81 

(2007). 

[17] G. Gulitski , M. Kholmyansky, W. Kinzelbach, B. Lthi, 
A. Tsinober and S. Yorish, J. Fluid Mech. 589, 83-102 
(2007). 

[18] G. Gulitski , M. Kholmyansky, W. Kinzelbach, B. Lthi, 
A. Tsinober and S. Yorish, J. Fluid Mech. 589, 103-123 
(2007). 

[19] L. Chevillard, C. Meneveau, L. Biferale, F. Toschi, Phys. 

Fluids 20, 101504 (2008). 
[20] P.E. Hamlington, J. Schumacher and W.J. A. Dahm 

Phys. of Fluids, 20, 111703 (2008). 
[21] L. Biferale, G. Boffetta, A. Celani, A. Lanotte, F. Toschi, 

Phys. of Fluids. 17(2), 021701/1-4 (2005). 
[22] D. Tordella and M. lovieno, J. Fluid Mech. 549, 441 

(2006). 

[23] D. Tordella, M. lovieno, S. Massaglia, Comput. Phys. 

Commun., 176, 539 (2007). 
[24] D. Tordella, M. lovieno, S. Massaglia, Comput. Physics 

Commun. 178, 883 (2008). 
[25] D. Tordella, M. lovieno, S. Massaglia, A. 

Mignone, in review for Comput. Phys. Commun., 

'http://arxiv.org/abs/1211. 1305^ (2013). 



